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ABSTRACT 

The long-term monitoring and high photometric precision of the Kepler satellite 
will provide a unique opportunity to sound the stellar cycles of many solar-type stars 
using asteroseismology. This can be achieved by studying periodic changes in the am- 
plitudes and frequencies of the oscillation modes observed in these stars. By comparing 
these measurements with conventional ground-based chromospheric activity indices, 
we can improve our understanding of the relationship between chromospheric changes 
and those taking place deep in the interior throughout the stellar activity cycle. In 
addition, asteroseismic measurements of the convection zone depth and differential 
rotation may help us determine whether stellar cycles are driven at the top or at the 
base of the convection zone. In this paper, we analyze the precision that will be pos- 
sible using Kepler to measure stellar cycles, convection zone depths, and differential 
rotation. Based on this analysis, we describe a strategy for selecting specific targets 
to be observed by the Kepler Asteroseismic Investigation for the full length of the 
mission, to optimize their suitability for probing stellar cycles in a wide variety of 
solar- type stars. 
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INTRODUCTION voke a combination of differential rotation, convection, and 

meridional flo w to modulate the field through a magnetic 

Despite the enormous efforts that have been carried out over ^ namo (e iDikpati fc Gilmanlbood ) 
the past 400 years to observe and understand the solar cycle, 

the delayed onset of cycle 24 suggests that w e still lack a Significant progress in dynamo modeling could only oc- 

comp lete understanding of this phenomenon (|Zimmermanl cur after helioseismology provided meaningfu l constraints 

120091 ). One of the best ways to improve our understanding ,° n Sun S inter i° r st ructure and dynamics (|Brown et all 

may be to study activity cycles in other stars, to help us Il989l: ISchou et al.lll998T). Generally the models assume that 

understand stellar cycles in general and not just the solar the driving of the solar dynamo is related to large gradients 

cycle. Asteroseismology is emerging as one of the best tools in the rotation rate of different layers of the Sun, presumably 

available for the study of stellar cycles. those layers with the largest gradients, i.e. a t the top and 

Astronomers have been making telescopic observations bottom of the near surface convection zone jBrandenburj 

of sunspots since the time of Galileo, gradually building a 120051). 

historical record showing a periodic rise and fall in the num- Stellar activity has been monitored in more than 100 

ber of sunspots every 11 years. We now know that sunspots stars ove r the last 40 years with the Mount Wilson sur- 

are regions with an enhanced local magnetic field, so this 11- vey JWilsonll978l : lBaliunas fc Vaughanll^ : lBaliunas^t~al] 

year cycle actually traces a variation in surface magnetism. 1 19951 ). This survey has revealed that around half of the 

Attempts to understand this behavior theoretically often in- solar-type stars show cl ear periodic cycles, w ith periods be- 
tween 2.5 and 25 years i|Baliunas et al.ll 19951 ). Later studies 
of this and other sampl es such as the Lowell Observatory 

* E-mail: y.p.elsworth@bham.ac.uk survey (|Hall et al,ll2007l ) suggest that there are two differ- 
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ent branches of stella r cycles in solar-type stars - one ac- 
tive and one inactive l|Saar fe Brandenburglll999l) . This bi- 
furc ation, which is also known from the Vaughan-Preston 
gap (jVaughan fc Prestonll 19801 ) between active stars younger 
than 2.5 Gyr and older less active stars, can also be seen in 
the lengths of the stellar cycles a s a function of the rota- 
tion period. IBohm- Vitensel (|2007h suggested, based on the 
fact that the number of rotations during a cycle seems to be 
different between the two subgroups of stars, that the bifur- 
cation is caused by different dynamos operating in the two 
branches. Active stars generally have around 300 rotations 
per cycle, whereas inac tive stars have around 100 rotations 
per cycle. According to IBohm- Vitensel (2007), this suggests 
that stellar cycles in active stars are generated by the gra- 
dient in the rotation rate close to the surface, whereas the 
stellar cycles in inactive stars are generated by the gradi- 
ent in the rotation rate at the base of the convection zone. 
In this scenario, a star with a given mass will start its life 
on the main sequence as a fast rotating active star with a 
relatively thin convection zone l|Bohm-Vitensell2007h . Since 
the convective turnover time in stars with thinner convec- 
tion zones is relatively short, these stars will have the largest 
rotational gradients near the surface, which would thus be 
where the dynamo would be driven in the stars on the ac- 
tive sequence. As the star evolves to the inactive sequence 
its surface rotation rate decreases significantly, which might 
be caused by additional deep m ixing in a growing convec- 
tion zone (|Bohm-Vitenseil2007h . This would result in the 
rotational gradient at the base of the convection zone grow- 
ing larger than the rotational gradient near the surface, and 
therefore the dynamo would be driven at the base of the 
convection zone in the inactive stars. The larger rotational 
gradient would also require fewer rotations for the dynamo 
to create a toroidal magnetic field strong enough to rise to 
the surface. 

Unfortunately, the Sun in many ways seems to fall be- 
tween the two sequences of active and inactive stars - es- 
pecially if one looks at the leng th of the stellar cyc le as a 
function of the rotation period ijBohm- Vitensel [ioorl l. This 
suggests that the Sun is not a good place to study stellar 
cycles, because two different dynamos might be operating 
in the Sun at the same time, confusing the picture. In other 
words, since the Sun seems to be in the transition from the 
active to the inactive sequence, it might be experiencing two 
different dynamos; one operating at the base of the convec- 
tion zone and one operating at the top - as opposed to most 
other solar-type stars which are only experiencing one of the 
two dynamos. Evidence of this scenario can also be found 
in the discussion of whether the solar dynam o operates at 
the b ase or at the top of the convection zone l|Brandenburg| 
l2005h . 

If we are to improve our understanding of the solar cy- 
cle, it seems that we need to study stellar cycles. To eval- 
uate the different dynamo models, we need asteroseismic 
measurements to complement the activity measurements of 
stellar cycles - both in the sense of asteroseismic measure- 
ments of cycle-induced changes in the oscillation mode fre- 
quencies and amplitudes, and in the sense of asteroseismic 
measurements of the convection zone depth and differential 
rotation. The problem is that the st ars in the Mount Wi l- 
son survey are relatively faint (V ~ 6: Bali unas et al1ll995h . 
making them unfit for dedicated ground-based asteroseis- 



mic networks like SONG (jGrundahl et al.l 120081 ). which are 
needed to obtain the extremely high frequency precision re- 
quired for the analysis. The reason for this is of course that 
bright G and K stars are rare in the sky, the few exceptions 
are: a Cen A & B (V = -0.01 and 1.33, respectively), 70 
Oph A (V = 4.03), k 1 Cet (V = 4.83) and the solar twin 18 
Sco (V = 5.5). It is therefore clear that the Kepler mission 
offers a unique possibility for sounding stellar cycles. 

Once we have obtained a firm understanding of the so- 
lar cycle, we might be able to use the fact that the Sun 
seems to have a more complicated dynamo than most other 
solar-type stars as an advantage - as this may enable us to 
test more complicated dynamo models on the Sun. But the 
prolonged delay in the onset of cycle 24 has clearly shown us 
that we do not hav e a firm understanding of the solar cycle 
l| Zimmerman! 120091 ) . 

This paper is the first in a series dedicated to the Ke- 
pler mission's study of stellar cycles in solar-type stars. The 
paper evaluates the expected frequency precision for Kepler, 
the changes in the behavior of the stellar cycles on the main 
sequence and the effect of the convection zone depth and 
differential rotation on the acoustic spectra. Although the 
analysis in this paper is dedicated to the Kepler mission, 
it can also serve as a guide for selecting targets for other 
asteroseismic surveys, such as SONG. 

We continue the analysis of the potential for soundin g 
activity cycles in solar- type stars bv lChaplin et all l|2007l ). 
but we now have the possibility to use specific estimates 
of the frequency precision for the Kepler mission. We also 
discuss problems with the way that t he acoustic amplitud e 
of the stellar cycles was estimated by IChapl in ct al. (2007). 
The analysis of the s ignature of stel l ar rota tion in the acous- 
tic spectrum follows IChaplin et al] (|2008aT ) , but here we ex- 
tend the analysis to include (1) updated estimates of mode 
lifetimes and the effect of these on the possibility of detect- 
ing the small rotational frequency splittings, and (2) the 
effect of differential rotation. Part of the analysis in this 
paper was done within asteroFLAG hare-and-hounds ex- 
ercises as part of th e preparation for the Kepler mission 
ijChaplin et al.ll2008bh and some of the new scaling relations 
presented in this paper have thus been obtained from and 
tested on the artificial stars generated as part of the Aster- 
oFLAG hare-and-hounds exercises. 

The layout of the rest of the paper is as follows. In 
section 2 we review the details of the study of solar-type 
stars under the Kepler Asteroseismic Investigation (KAI). 
A review of the planned ground-based support for the study 
of stellar cycles in solar-type stars under the KAI is given in 
section 3 and an analysis of the expected frequency precision 
is given in section 4. Sections 5-7 use a number of isochrone 
tracks to estimate the effect of rotation (section 5), stellar 
cycles (section 6) and the convection zone (section 7) on 
the oscillation mode frequencies that we expect to measure 
from the observations. This analysis is used in a discussion 
of the best strategy for selecting the optimal full-length of 
the mission Kepler targets to increase our understanding of 
stellar cycles (section 8). Concluding remarks are found in 
section 9. 
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Figure 1. Padova isochrone tracks used in the analysis, shown 
as a function of effective temperature and surface gravity. The 
different line styles give the age of the stars: the solid line repre- 
sents stars that are f .0 Gyr old, dotted f .6 Gyr, dashed 2.5 Gyr, 
dot-dashed 4.0 Gyr, dot-dot-dot-dashed 6.3 Gyr and long dashed 
fO.O Gyr. This notation is used for the rest of the paper. The 
gray dots are stars from the Kepler Input Catalog with g-band 
magnitudes between 8 and fO. 



2 KEPLER OBSERVATIONS 

Kepler is a NASA Discovery class space mission that was 
launched into an Earth trailing orbit on 7 March 2009, 
with the primary g oal of detecting te rrestrial planets using 
the transit method (Basri ct al. 200§|). The planned mission 
lifetime of Kepler is 42 months (excluding commissioning) 
with an optional extension of 30 months. The satellite is ex- 
pected to be operational beyond the extended mission, but 
the downlink capability will be limited due to the distance 
between the satellite and the Earth. Kepler will monitor over 
100,000 stars with a cadence of 30 min. In parallel with the 
planet pro gram, Kepler also has an Asteroseis mic Investiga- 
tion (KAI; lchristensen-Dalsgaard et alj|2007h . 

The KAI will begin with an initial run of 8 - 10 months 
depending on the length of the commission phase. In this 
survey phase of the mission, targets will only be observed for 
one month each. In the survey phase, close to 40,000 pixels 
are available for the KAI in short cadence (60 s). After the 
survey phase the KAI will shift to a specific target phase for 
the rest of the mission, and it is among these targets that it 
will be possible to study stellar cycles. In the specific target 
phase, the KAI will have the possibility to allocate 11,800 
pixels for observations in short cadence. A conservative es- 
timate is that around 3,000 of these 11,800 pixels will be 
allocated to F, G, and K stars on the main sequence suit- 
able for studying stellar cycles. The reason for giving these 
numbers in pixels and not as a quantity of targets is that 
bright stars will require more pixels than faint stars due to 
saturation. It is estimated that a 6th magnitude star will re- 
quire 2,640 pixels, whereas a 12th magnitude star will only 
need 70 pixels. In practice, this means that the specific tar- 
gets for studying stellar cycles could be chosen as e.g.: one 
6th magnitude, ten 8th magnitude or 43 13th magnitude 
stars. The numbers of pixels that need to be allocated for 
stars at a given brightness are likely to decrease after the 
survey phase, when the characteristics of the instrument are 
better known. Kepler is designed to provide nearly photon 



Figure 2. All stars in the Kepler Input Catalog with g-band 
magnitudes between 8 and 10 (gray dots) and targets chosen for 
observations in the first three months of the survey phase (black 
crosses). 

noise limited photometry on saturated stars. The concept of 
doing photometry on sat urated stars ha s been demonstrated 
with HST observations l|Gillilandl200sh . 

The targets for the su rvey phase have bee n chosen from 
the Kepler Input Catalog (jLatham et al.ll2005l ). and specific 
targets to be observed for the full-length of the mission must 
be chosen from the targets in the survey. It is expected that 
only observations from the first three months with be avail- 
able before the specific targets need to be chosen, so in prac- 
tice only targets from the first three months of observations 
can be selected as specific full-length of the mission targets. 
Fig. 1 shows all of the stars in the Kepler Input Catalog with 
g-band magnitude between 8 and 10, with Padova isoc hrone 
tracks i|Bonatto et alj|2004 [Girardi et al.ll2002l . \20o4 ) plot- 
ted on top as a function of effective temperature and surface 
gravity. The isochrone tracks were calculated for 6 different 
ages between 1 and 10 Gyr in steps of 0.2 dex, using a metal- 
licity of z = 0.02. The same stars are plotted in Fig. 2, but 
here the targets selected for observations in the first three 
months of the survey phase are marked. It is seen that F 
stars are highly over-represented compared to cooler stars 
in the first three months of the survey phase and in general. 
Although this reflects the stellar populations at a given mag- 
nitude bin, it is unfortunate as F stars are the ones that are 
hardest to analyze with asteroseismolog y since they d o not 
obey the asymptotic freq uency relation ijTas soul 1980) and 
have short mode lifetimes (|Chaplin et al.|[2009l ) . The fraction 
of stars which show stellar cycles is also sig nificantly larger 
amon g G and K stars than among F stars (Ba liunas et al.l 
ll995l V 



3 GROUND-BASED OBSERVATIONS 

Complementary ground-based chromospheric observations 
are important for mainly two reasons: firstly, they can be 
used to compare the cycle-induced changes in the chromo- 
sphere to changes in the amplitudes and frequencies of the 
oscillation modes - in other words, the complementary data 
provide a possibility to link the changes near the surfaces of 
the stars (i.e. in the chromospheres) to changes in the stellar 



4 C. Karoff et al. 



10 r 




0.1 1 F 



6500 6000 5500 

T.„ [K] 



5000 



4500 




06 ^ ......... , 

6500 6000 5500 

T.„ [K] 



5000 



4500 



Figure 3. Mode height as a function of effective temperature 
for the Padova isochrone tracks. It is seen that the mode height 
increases dramatically when the stars begin to evolve off the main 
sequence. 



interiors (i.e. in the oscillation modes). Secondly, ground- 
based chromospheric observations will be important if the 
cycles are significantly longer than the Kepler mission life- 
time, or if the cycle-induced changes in the mode frequencies 
and amplitudes are too small to be measured with seismol- 
ogy. We will analyze the expected cycle lengths and seismic 
amplitudes in section 6. 

We plan to measure the chromospheric activity of the 
coolest candidate targets discussed in section 2 during the 
summer of 2009 using the high-resolution Fibre-fed Echelle 
Spectrograph ( FIES) mounted on the 2.6 m eter Nordic Op- 
tical Telescope (|Frandsen fc Lindberd l^OOO - ), After the spe- 
cific targets have been selected, we will start to monitor the 
targets to look for stellar cycles in those stars. Note that 
even stars which show no sign of stellar cycles will be in- 
teresting, since it might be possible to use the asteroseismic 
measurement of the depths of the convection zones and ro- 
tation to understand why these stars show no sign of stellar 
cycles. 



4 EXPECTED FREQUENCY PRECISION 

To estimate the potential for sounding stellar cycles, and for 
measuring the depths of the convection zones and differential 
rotation in the stars observed by Kepler, we first need to 
know what frequency precision we will be able to obtain from 
the observations. The frequency precision is determined by 
three parameters: the S/N of the oscillation modes in the 
acoustic spectrum, the mode lifetime r and the length of 
the ob servations T as given by the formulation of iLibbrechtl 
<|1992T ): 



A-kT ' 



(1) 



where F(f3) is a function of the inverse S/N, /3: 

F(/3) = y/(l + 0) + + \//?) 3 , (2) 

and A = 1/(ttt) is the mode FWHM linewidth. The Lib- 
brecht formulation provides a lower limit on the frequency 
precision, since it does not account for blending effects in 



Figure 4. Frequency precision after 1 year of observations of a 
9th magnitude star with Kepler as a function of effective temper- 
ature for the Padova isochrone tracks. The frequency precision 
scales with l/\/T where T is the length of the observations. The 
frequency precision after 4 years of observations can therefore be 
calculated by dividing the numbers on the y-axis by 2. 



closely spaced oscillation modes. This is mainly a problem 
for modes with different azimuthal order, which is analyzed 

in detail in section 5. 

W e use the T e g 4 scaling relation from IChapl in et al] 
(2009) for the mode lifetimes, which also provides us with 
the following relation for the height H of the oscillation 
modes in the acoustic power density s pectrum by using the 
scalin g relation for the amplitude from lKieldsen fc Beddina 
|l995l ) for intensity observations: 



H oc g~ 



(3) 



where g is the surface gravity. Using the Padova isochrone 
tracks, we have plotted the mode height as a function of 
effective temperature in Fig. 3. It is seen that the mode 
height increases dramatically when the stars begin to evolve 
off the main sequence. It is also seen that the mode height 
decreases toward lower temperature, from 3 ppm 2 //iHz for 
the Sun down to just over 1 ppm 2 //iHz for stars with tem- 
peratures around 4500 K. 

Knowing the expected mode heights, we can use the ex- 
pected photometric errors of the Kepler photometer to cal- 
culate the expected S/N. The photometric error is a function 
of many factors - the main one being the magnitude. It is 
estimated that the photometric error for a 60 sec integration 
(Tint with Kepler will be 45 ppm at 8th magnitude, 69 ppm 
at 9th magnitude and 111 ppm for a 10th magnitude star. 
The photometric errors are based on laboratory tests with 
fully flight-like electronics with an assumed lev el of stellar 
(activ ity and granulation) noise of 10 ppm (see lKoch et al] 
|2000| . for details of the laboratory tests) . Assuming that the 
noise is random and uncorrelated, the photometric error can 
be converted into a power density error <j p d using: 



2Ti nt a in 



(4) 



where T ln t is the sampling time - i.e. 60 sec. In this way, we 
expect the error in the power density spectrum to be 0.24, 
0.57 and 1.48 ppm 2 //dTz for a 8th, 9th and 10th magnitude 
star, respectively. In comparison MOST, WIRE and CoRoT 
have obtained noise levels in the high-frequency part of the 
power density spectrum near the Nyquist frequency as low 
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Figure 5. Rotation frequency splitting as a function of effective 
temperature for the Padova isochrone tracks. It is seen that the 
rotational splitting decreases toward lower effective temperature. 



as 3. 1, and 0.15 ppm 2 / / iHz, respectively jGuenther et all 
|200S| ; iBruntt et ail 120051 ; lAppourchaux et alj|2008h . These 
values are not directly comparable to the expected error in 
the power density spectrum that we calculate for Kepler 
since they are measured at high frequency and thus include 
very little activity and granulation signal, but it is clear 
that the expected error that we calculate does not seem 
unrealistic in comparison. 

The estimated frequency precision for a 9th magnitude 
star as a function of effective temperature is shown in Fig. 4. 
It is seen that the evolved stars will yield the best frequency 
precision. It is also seen that the frequency precision im- 
proves toward lower effective temperature, from around the 
effective temperature of the Sun down to an effective tem- 
perature of 4500 K, but the slope is small. In other words, 
although the mode height changes on the main sequence, 
the frequency precision for oscillation modes in stars on the 
main sequence is approximately constant. 



Figure 6. Ratio between rotational frequency splitting and the 
mode linewidth as a function of effective temperature for the 
Padova isochrone tracks. Note the interesting conclusion that the 
effective temperature of the Sun is the worst possible effective 
temperature at which rotational frequency splittings can be mea- 
sured. This provides some hope that we will be able to measure 
differential rotation with asteroseismology. 



As the oscillations that we are studying in this paper are 
high-order p modes we can also assume that the rotation 
frequency splittings equal the rotation rates. Note that this 
is not the case for, e.g. the high-order g modes in white 
dwarfs where th e rotational splitting is only equal to half 
the rotation rate l)Christensen-Dais gaard 2003]). Under these 
assumptions the rotation period can easily be converted into 
an equivalent rotational frequency splitting by using the re- 
lation: Prot = 1/Afrot- 

Fig. 5 shows the rotational frequency splitting as a func- 
tion of effective temperature for the Padova isochrone tracks. 
It is generally seen that the rotational splitting deceases to- 
ward lower effective temperature and that the slope of this 
trend is much larger for stars hotter than the Sun compared 
to stars cooler than the Sun. 



5 PREDICTION OF ROTATION 

The mean surface rotation period P Iot can be estimat e d from 
the empirical relationship given by lAigrain et al.l l|2004l ). 
This relationship is de rived from photometric o bservations 
of stars in the Hyades l|Radick et al.lll987l . [l995T ). The rela- 
tionship is: 

log(P rot ) - 0.51og(i/0.625) = 

f -0.669 + 2.580(5 - V) for 0.45 < B - V < 0.62 (5) 
\ 0.725 + 0.326(B - V) for 0.62 < B - V < 1.30 

Here, t is the age of the star in Gyr, and it ap pears in the 
relatio nship because the t 1 ^ 2 spin-down law of ISkumanichl 
1 19721 ') is implicit. 

When measuring stellar rotation with asteroseismology, 
it is the rotational frequency splitting Az/ rot that is mea- 
sured rather than the rotation period. Assuming that the 
stars possess internal rotation rates that are comparable 
to the surface rates, as in the Sun, it is reasonable to as- 
sume that the mean rotation rate is given as the inverse 
of the mean rotation period. This might be questionable 
for younger stars fsee lChaplin et al]|2008al . for discussion). 



5.1 Effect of finite mode lifetime 

The possibility of measuring the rotational frequency split- 
ting, and thus the mean rotation period of the stars, relies 
not only on a frequency precision sufficient to measure the 
small splittings, but also on splittings that are relatively 
large compared to the mode linewidth. If the linewidth is 
significantly larger than the splitting, it will not be possi- 
ble to disentangle the modes of different azimuthal order 
m, and it will thus not be possible to measure the splittings. 
For this reason the internal rotation rate of the Sun was first 
obtained from low-degree helioseismology when it was real- 
ized that only low-frequency m odes with long mode lifetime 
should be used for this analysis l|Elsworth et al.|[l995l ). Fig. 6 
shows the ratio between the rotational frequency splittings 
and the mode linewidth, which can serve as a measure of the 
chances of disentangling the modes with different azimuthal 
order. The figure shows the surprising result that the effec- 
tive temperature of the Sun is the worst possible effective 
temperature at which rotational frequency splittings can be 
measured. 



6 C. Karoff et al. 



As shown by iGarcfa et alj (|2004l ). it is just possible 
to measure radial differential rotation in the low-order fre- 
quency splittings from 6 years of low-degree observations of 
the Sun. Since the ratio between the rotational frequency 
splitting and the mode linewidth will be larger for all other 
stars on the main sequence, it is therefore not unreasonable 
to hope that we will be able to use Kepler to measure radial 
differential rotation in some solar-type stars with rotation 
profiles similar to the Sun. 

In order to measure stellar rotation the rotation axis 
of the star needs to be inclined with respect to the line-of- 
sight. If the inclination angle between the rotation axis of the 
star and the line-of-sight is low, then the rotation will not 
significantly affect the frequencies of the observed oscillation 
modes with different azimuthal order, and the rotational fre- 
quenc y splitting can thus not be measured. iGizon fc SolankH 
( 2003) find that it is only possible to measure the rotational 
frequency splitting if the inclination angle is larger than 
around 30°. Statistically, assuming random orientation of 
the rotation axis, it will therefore not be possible to mea- 
sure the rotational frequency splitting in around 15 per cent 
of the targets due to low inclination angles. 

5.2 Effect of differential rotation 

The Sun rotates fastest in the outer part of the convection 
zone beneath the equator, with a rotation rate of around 
470 nHz. The rotation rate decreases to around 435 nHz at 
the base of the convection zone, and to around 450 nHz at 
the surface - with a further decrease along the surface to the 
poles with a value of around 370 nHz at 60° which is the 
highest latitude at which rot ation rates can be measured re- 
liably with helioseismology (jSchou et al.lll998h . The largest 
gradients in the rotation rate are generally seen near the 
surface and near the base of the convection zone. 

Differential rotation in the Sun and in other solar-type 
stars will manifest itself both as latitudinal and radial dif- 
ferential rotation. Latitudinal differential rotation is mainly 
on the surface, and can be measured either with asteroseis- 
mology or by fitting a spot model to the observed intensity 
time series as has been don e for the stars k 1 Ceti and e Eri 
using the MOST satellite (|Croll et al.l 120061 : IWalker et all 
120071 ). Radial differential rotation is confined to the stellar 
interior, and can only be measured with asteroseismology. 

Differential rotation can be measured with asteroseis- 
mology by measuring differences in the separation between 
modes with different azimuthal order m (the rotation split- 
ting) over different radial orders n and angular degrees I. 
A simplified way of doing this is by averaging the rota- 
tional splittings over different radial orders, and thus ob- 
tainin g the splitting only a s a function of the angula r degree 
(see lGizon fc Solankilliool . 120041 : iBailot et al.ll200fjh . This is 
done to increase the S/N of the measured splittings, but the 
approach in strongly biased towards measuring latitudinal 
differential rotation on the surface, and is much less efficient 
at measuring radial differential rotation at the base of the 
convection zone. For the Sun, the mean difference in the ro- 
tational splitting betwee n 1 = 1 and 1 = 2 m odes is around 
7 nHz at low frequency {Garcia et alj [20031 ). We therefore 
adopt this as a conservative reference value, although there 
might be better ways of averaging the rotational splittings 
for stars with thicker convection zones and thus larger radial 
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Figure 7. The difference in the rotational splitting between 1 = 1 
and I = 2 modes averaged over different radial orders as a function 
of effective temperature for the Padova isochrone tracks. It is 
seen that the differences increase toward higher temperature and 
decrease with age. 

rotation gradients (|Bohm- Vitensdl2007h . Such new averag- 
ing schemes would need to be calculated based on detailed 
forward modeling of the expected changes in the rotational 
splitting, not only over different angular degrees, but also 
over different radial orders. In this way, it would hopefully be 
possible to measure average rotational splittings that would 
be more sensitive to radial changes in the rotation rate at 
the base of the convection zone than to latitudinal changes 
in the rotation rate at the surface. 

Since only surface (latitudinal) differential rotation has 
been measured in other solar-type stars so far, we are only 
able to make predictions of the changes in the surface dif- 
ferential rotation on the main sequence. Surface differential 
rotation is normally expressed as the change in the rota- 
tion period AP ro t over the stella r disc. A relation betw een 
AP r ot and P rot was obtained by iDonahue et al.l (|l996T ) by 
looking at the changes in the rotation period measured in 
at least five different observing seasons in the stellar Ca II 
measurements from Mount Wilson Observatory: 

AP rot oc P^ ±0 '\ (6) 

The interesting parameter for asteroseismology is not the 
change in the rotation period AP rot caused by differential 
rotation, but rather the change in the rotation rate AQ. 
This can easily be obtained by assuming that AP ro t <C Prot 
and thus Afi oc AP ro t/P r 2 ot , which gives the following re- 
lation between the changes in the rotation rate caused by 
differential rotation and the mean rotation period: 



Adopting 7 nHz as the mean difference in the rotational 
splitting between I = 1 and I = 2 modes at low frequency 
{Garcia et alj |2003) in the Sun, we can plot the mean differ- 
ence shown in Fig. 7. It is seen that the differences increase 
toward higher temperature and decrease with age. 

The result in Fig. 7 should be seen as a lower limit on 
the effects of differential rotation that can be detected with 
asteroseismology, both because there might be better ways 
to measure radial differential rotation (as discussed above) 
and because differential rotation measured through changes 
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Figure 8. Stellar cycle period as a function of effective temper- 
ature for the Padova isochrone tracks. No general trend is seen. 



in t he rotation period sa mples a limited range of latitudes 
fsee lDonahue et al.lll996l . for discussion). 

Slightly hi gher values o f the e xponent in Eq. 6 have 
found by Henry et all (|l995h . who obtained an ex- 
ponent of 1.76 ± 0.06 based o n photometry of active bi- 
naries, and I Barnes et ahl (|2005l ) who obtained 1.85 ± 0.10 
based on Doppler imagin g techniques and spectral analysis 
(R einers fc Schmitt 2003). We have chosen to use the rela- 
tion by iDonahue et alj l|l996i) because the sample used in 
our study is not biased towards really active stars (as the 
photometry results) or towards really fast rotating stars (as 
the Doppler imaging and the spectral analysis results). The 
effect of increasing the exponent in Eq. 6 would be to flatten 
the curves in Fig. 7. 

Differential rotation in solar- type stars has also bee n 
studied with mean- field models bv lKiiker fc Riidigerl |2008), 
who see a similar relation between Af2 and the effective 
temperature as we see in Fig. 7 - i.e. the slope of the curves 
seems to increase dramatically at an effective temperature 
near 6000 K. 



6 PREDICTION OF STELLAR CYCLES 

A commonly used indicator of surface activity on stars is 
the Ca II H and K emission index. This index is usually 
expressed as i? H K ~ the average fraction of the stellar lumi- 
nosity tha t is emitted in t he Ca II H and K emission line 
cores (see lHall et aL 1 120071 , for dis cussion). Using the first 
data from the Mount Wilson survey, iNoves et all jl984h wort- 
able to obtain relations between the rotation period P ro t , the 
R'hk index and the B — V color, which can be used to obtain 
a scaling relation for the Phk index: 

log(r c /P rot ) = -(0.324 - 0.400y + 0.283y 2 - 1.325y 3 ), (8) 

where y =log(P H K x 10 ) and r c is the convective turnover 
time. The convective turnover time may then be calculated 
from: 



log (t c ) = 

1.361 + 0.166a; +0.025:r 2 
1.361 - 0.140a; 



5.323^ 



for 
for 



x ^ (9) 
x < 



Figure 9. Stellar cycle frequency shifts as a function of the 
effective temperature for the Padova isochrone tracks. For the 
black lines the shift s were calculated usi ng the scaling relation 
8u cyc oc AR' nK from IChapl in et al. I J2007t> and for the gray lines 
the shifts were calculated using the scaling relation 8u cyc oc 
fl 2. 5^0.26 Afl ^ from iMetcalfe et al.l d2007h . It is seen that the 
shifts calculated using the scaling relation by IChaplin et al.l 
i2007f ) are nearly constant for stars cooler than 5500 K, while 
at the same age the shi fts calculated using the scaling relation 
by IMetcalfe et al] J2007T) generally decrease around 50 per cent 
between 5500 K and 4500 K. 



cubic Eq. 8 using the known input parameters: B — V and 

Prot- 

Knowing the -Rhk index, we ca n calculate the c ycle pe - 
riod P cyc from the relation given bv lChaplin et al l j2008a): 



log(Pcyc) = 

-6.7 - 1.61og(P HK ) 
-5.5 - 1.41og(P HK ) 



for log(P HK ) < -4.75 (10) 
for log(P HK ) > -4.75 



Fig. 8 shows the expected stellar cycle period as a func- 
tion of effective temperature. No general trend is seen. For 
the three youngest isochrones, it is seen that the cycle pe- 
riod generally increases toward higher temperature. For ef- 
fective temperatures lower than the effective temperature of 
the Sun, it is seen that the cycle periods generally increase 
with age as the stars evolve from active to inactive stars 
while their rotation periods decrease. On the other hand, 
as the stars evolve from the active to the inactive sequence, 
the number of rotations in each stellar cycle decreases and 
therefore no general trend is seen in the figure. 

The amplitudes of the activity cycles AR' HK can be 
obtained from the relation given by Saar & Brandenburg 
(2002): 



APhk oc (P H k)° 



(11) 



where x = 1 — (B — V) and the convective turnover time is 
measured in days. To calculate the Phk index, we solve the 



The p-mode frequencies in the Sun are generally believed to 
be affected by the solar cycle through changes in the tur- 
bulent velocity in the outermost regions of the near surface 
convection zone, caused by changes in the RMS m agnetic 
field in the Sun (|Dziembowski fc Goode|[20o3 , 120051 ). Rela- 
tions between the Ca II amplitudes of the activity cycles 
A.Rhk a nd the cycle fr e quenc y shi fts Sv cvc have been p ro- 
pos ed bvlChaplin et al] d2007h and IMetcalfe et al] (|2007t >. 

IChaplin et al] 1 2007f ) proposed that the cycle frequency 
shifts simply scale linearly with the amplitudes of the activ- 
ity cycles: 
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Su cyc OC Ai?H 



(12) 



Evidence for this relation is found from observations of low- 
degree oscillations in the Sun, where it is seen that the cycle 
frequency shifts change linearly with chang es in the Mg II 
lines over the solar cycle (Cha plin et al.ll2007l ). The expected 
stellar cycle frequency shifts as a function of the effective 
tempe rature obtained using the scaling from IChapl in et al.l 
is shown in Fig. 9 (black lines). It is seen that the 
cycle frequency shifts generally increase toward lower tem- 
perature and decrease with age. 

On the other hand, iMetcalfe et all (|2007f) assume that 
the cycle frequency shifts scale as 



Su cvc oc ^A_Rt 



(13) 



where D is the depth of the source of the perturbations be- 
neath the photosphere and I is the mode inertia. The inverse 
relation with the mode inertia is kn own from studies of so- 
lar cycle mode frequency shifts (se e Libbrecht fe WoodardI 
199d; iDziembowski fc Goodel 12004 bo05l . for discussion). 
Metcalfe et alj (|2007I ) then assume that the mode inertia 
scales as M/R and that the source depth scales with the 
pressure scale height at the photosphere H p , which itself 
scales as L°- 25 R 15 /M. We thus obtain: 



n2.5 rO.25 
O^cyc OC — Ai? H 



(14) 



Note that we lack a physical argument for the relation be- 
tween the source depth and the pressure scale height other 
than the fact that most things in the near surface layers 
scale with the pressure scale height. The Kepler observa- 
tions might eventually provide us with such arguments as it 
might be possible to extract information of the location of 
the source of the p-mode excitation from the high-frequency 
oscill ations in the solar- type stars as has been done for the 
Sun iKumar fc Ltj||l99lh and other solar-type stars (|Karofl1 
12007ft , but it remains to be proven if similar kinds of analysis 
could also reveal the location of the source of the perturba- 
tions caused by the stellar cycles. 

Fig. 9 also shows the expected stellar cycle frequency 
shifts as a function of the effective temperat ure obtained by 
using the scaling from IMetcalfe et al.l (|2007f ) (gray lines). It 
is seen that the shifts generally increase toward higher tem- 
perature and decrease with age. The shifts calculated using 
the scaling relation by I Chaplin et all (|2007f ) are nearly con- 
stant for stars cooler than 5500 K but with the same age. 
On the ot her hand, the shifts c alculated using the scaling 
relation by IMetcalfe et ail l|2007l ) generally decrease around 
50 per cent between 5500 K and 4500 K for a given age. 
It is currently not possible to determine which of the two 
scaling relations is more appropriate, but it is clear from 
Figs. 4 & 9 that we will be able to measure the oscillation 
mode frequencies precisely enough from one year of Kepler 
data to detect stellar cycle frequency shifts when these fre- 
quencies are compared to Kepler observations obtained in 
other years. We also note that the Kepler observations might 
allow us to test which of the scaling relations is the better 
one. 



7 PREDICTION OF THE CONVECTION 
ZONE DEPTH 

One of the most important things that will help us under- 
stand the stellar cycles we observe with Kepler is a mea- 
surement of the convection zone depth. The depth of the 
convection zone can be measured with asteroseismology by 
analyzing the second differences in the frequencies: 



AoIAi 



n + l,l 



(15) 



Localized regions of rapid variation of the sound speed, 
such as the bottom of the convection zone, will induce 
an oscillatory component in the se cond frequency differ- 
ences as a function o f frequency (|Monteiro et al.l l200d : 
iHoudek fc Goughll2007r i. Other localized regions that cause 
oscillatory components in the second fre quency differences 
are th e ionization zones of He I and He II. IHoudek fc Goughl 
have given full analytical solutions for how to rep- 
resent the oscillatory components in the second frequency 
differences from the three regions. For the analysis here, we 
are only interested in determining how the amplitude of the 
oscillatory component originating from the base of the con- 
vection zone depends on fundamental stellar parameters. 

To evaluate how the amplitude of the oscillatory com- 
ponent originating from the base of the convection zone 
changes with fundamental stellar parameters, we have 
used frequencies from 6 artificial solar-type stars gener- 
ated as part of the as teroFLAG hare-and-hounds exercises 
|Chaplin et all l2008bh . The amplitudes of the oscillatory 
components in these stars were measured by calculating the 
least-squares periodogram of the second differences. These 
periodograms all showed 3 peaks from the 3 regions of rapid 
variation of the sound speed, and the peak with the high- 
est frequency is the one that originates from the base of the 
convection zone. We then looked for correlations between 
the amplitude of the peak originating from the base of the 
convection zone and the effective temperature and surface 
gravity. It was seen that there was a clear correlation be- 
tween the logarithm of the surface gravity log(gr) and the 
amplitude of the form: 



A cz oc log (g). 



(16) 



We adopt a value of 0.4 ^tHz for the mean amplitude of the 
oscillator y component in the seco nd frequency differences in 
the Sun (|Houdek fc Goughl l2007h . Fig. 10 shows the ampli- 
tude of the oscillatory component originating from the base 
of the convection zone as a function of the effective tem- 
perature. It is seen that the amplitude generally increases 
tow ard lower t e mpera ture. 

iRoxburghl {2009) has recently suggested another ap- 
proach for measuring the convection zone depth with aster- 
oseismology. This method makes use of the second and third 
frequency differences between modes of angular degree and 
1 instead of the second differences of modes with the same 
degree. This method has the advantage that the amplitude 
of the oscillatory component of the differences is constant as 
a function of frequency - whereas it decreases toward higher 
frequencies for the differences from modes with the same de- 
gree. Also, the method can be tuned so that the signal from 
the He II ionization zone cancels out, which makes it easier 
to detect the signal from the base of the convection zone. 
However, the method has two weaknesses compared to the 
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Figure 10. The amplitude of the oscillatory component in the 
second frequency differences originating from the base of the con- 
vection zone as a function of the effective temperature for the 
Padova isochrone tracks. 



differences from modes with the same degree. Firstly, it only 
uses information from modes with degree and 1, and not 
from modes with degree 2 that we will also detect in the 
Kepler observations (which means that not all information 
from the observations is used to measure the depth of the 
convection zone). Secondly, it is not clear whether the signal 
from the He II ionization zone can be canceled out in stars 
other than the Sun. 

The si gnal from the bas e of the convection zone in the 
method by Roxburgh! l|2009l ) will also scale with log(g) as 
suggested in Eq. 16, though the absolute amplitude and S/N 
could be different from the amplitudes and S/N obtained for 
the second differences from modes with the same degree. 



hotter and cooler stars will have larger rotational splittings 
relative to the mode linewidths in the acoustic spectrum. 
Differential rotation will generally be measured most easily 
in young stars. For stars hotter than 6000 K, differential 
rotation will be measured most easily in the hottest stars, 
while there is almost no dependence on temperature for stars 
cooler than 6000 K. 



8.3 Optimal stellar parameters for activity cycles 

The analysis did not provide a clear answer to where on the 
main sequence it will be easiest to sound stellar cycles. It 
is clear that young stars will generally have shorter periods 
and larger frequency changes caused by stellar cycles. In 
order to use ast eroseismology t o inve stigate the hypothesis 
put forward by iBohm-Vitensel l|2007l ) - that the Vaughan- 
Preston gap is caused by two different kinds of dynamos 
operating on either side of the gap — it is important that 
stars on either side of the gap are observed with Kepler. The 
best way to determine which side of the gap the stars are on 
is by determining their ages (see section 8.6). 



8.4 Optimal stellar parameters for measuring the 
convection zone depth 

It will generally be easier to measure the depth of the stellar 
convection zone in cooler stars, since the amplitudes in the 
second frequency differences are expected to be larger than 
in hotter stars. The analysis has shown that we can expect to 
measure the convection zone depth in all KAI full length of 
the mission specific targets in the analyzed range in effective 
temperature. 



8 DISCUSSION 

Below we summarize the different recommendations for se- 
lecting asteroseismic targets for sounding stellar cycles. 

8.1 Effect of stellar brightness and S/N of 
observations 

At 8th magnitude the S/N of the oscillation modes will al- 
low the detection of modes in all stars, and the frequency 
precision will mainly be limited by (1) the mode lifetime, 
which will cause the power to be smeared out over a larger 
frequency range, and (2) the length of the observations. 9th 
magnitude stars will allow the detection of oscillations in 
most cases. Here the frequency precision for stars with long 
mode lifetimes will be limited not only by the mode lifetime 
and the length of the observations, but also by the S/N. 10th 
magnitude stars will only allow the detection of oscillation 
modes for temperatures similar to the Sun and hotter, and 
the frequency precision will be limited by the S/N for most 
mode lifetimes. 

8.2 Optimal stellar parameters for measuring 
stellar rotation with asteroseismology 

The worst place in parameter space for measuring stellar 
rotation with asteroseismology is the place of the Sun - both 



8.5 General optimal stellar parameters for 
sounding stellar cycles 

From these recommendations it is generally seen that the 
targets should be as young as possible, with the reserva- 
tion that both sides of the Vaughan-Preston gap need to be 
sampled. 

The analysis of the effect of differential rotation and 
stellar cycles on the acoustic spectra favors hotter stars, 
though the dependence is not strong and unambiguous. The 
effect of the convection zone depth favors cooler stars. We 
thus conclude that stars both cooler and hotter than the 
Sun should be on the KAI full length of the mission specific 
targets list, and to ensure a relatively equal distribution this 
means that we should favor cooler stars since they are clearly 
under- abundant among the targets already selected for ob- 
servations in the first three months of the survey phase. 

It is clear that measuring stellar differential rotation will 
be the most demanding task of those analyzed here, and this 
task will set high demands on the frequency precision. This 
means that it will only be possible to measure stellar differ- 
ential rotation with asteroseismology on Kepler stars after 
they have been observed for the full length of the mission. 
On the other hand, the analysis has shown that evaluating 
the potential for measuring stellar differential rotation based 
on the frequency shifts seen in the Sun from differential ro- 
tation is probably too pessimistic. This provides some hope 
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that we will be able to measure stellar differential rotation 
in some of the stars observed with Kepler. 

The analysis has also shown that it is probably going 
to be possible to measure the oscillation mode frequencies 
precisely enough in one year of Kepler data to be able to 
see stellar cycle frequency shifts when these frequencies are 
compared to Kepler observations obtained in other years. It 
also seems possible to measure the depth of the convection 
zone from just one year of Kepler observations, though the 
full data set will eventually be available for this purpose. 

It is also worth noting that mode frequency shifts might 
not be the only way to sound the stellar cycles. The stellar 
cycles will also manifest themselves a s changes in the am - 
plitudes and lifetimes of the modes (|Chaplin et al.| [200C0, 
and possibly also in flare induced high-frequency oscillations 
jKaroff fe Kieldsenll2008h . though the last concept remains 
to be proven. 



8.6 Asteroseismic determination of stellar age 

Since one of the recommendations is that mainly young 
stars should be selected, we need a way to determine the 
stellar ages. The best tool for doing this is asteroseismol- 
ogy, and it is reasonable to assume that a rough age esti- 
mate can be obtained from the survey observations by plot- 
ting the small and large frequency separations in the as- 
teroseismic HR-diagram (also kno wn as the JC-D diagram, 
IChristensen-Dalsgaard et al.l [2007) . Since we are interested 
in stars around the Vaughan-Preston gap (age 2-3 Gyr), we 
can safely assume that the lines in the asteroseismic HR- 
diagram are vertical to within the precision needed for Ke- 
pler specific target selection. This basically means that we 
are mainly interested in stars which have a small separation 
equal to the small separation in the Sun (6 /iHz) and larger. 



8.7 Recommendations for Kepler target selection 
for sounding stellar cycles 

We finally end up with the following selection procedure: 

(1) Use the targets with magnitudes between 8 and 10 in 
the first 3 months of the Kepler survey phase as candidates. 

(2) Start selecting stars from the cool end. 

(3) Ensure that the oscillations can be seen in the acoustic 
spectrum. 

(4) Ensure that the oscillation modes can be understood 
in the framework of the asymptotic frequency relation. 

(5) Ensure that hints of rotational splitting can be seen. 

(6) Ensure that the small separation is large (6 /iHz). 

(7) Ensure that both active and inactive stars have been 
selected (based on ground-based observations of chromo- 
spheric activity in the stars). 



9 CONCLUSIONS 

We have presented an analysis of the potential for using 
observations from the Kepler Asteroseismic Investigation to 
improve our understanding of solar and stellar cycles. The 
analysis has shown that it will probably be possible to sound 



stellar cycles with Kepler in stars brighter than 9th magni- 
tude by measuring the changes in the oscillation mode fre- 
quencies caused by the stellar cycles. It should be possible 
to measure the depth of the convection zone and the mean 
rotation rate from Kepler observations of stars brighter than 
9th magnitude. The analysis has also revealed that measur- 
ing stellar differential rotation will be highly demanding on 
the frequency precision, and it is clear that it will only be 
possible after the mission has finished. On the other hand, 
one of the unexpected results of the analysis was that the 
worst place in parameter space for measuring stellar rota- 
tion with asteroseismology is the place of the Sun - which 
provides some hope that it might be possible to measure stel- 
lar differential rotation in some of the stars observed with 
Kepler. 

The scaling relations we use for rotation, mode lifetime, 
activity etc. represent a statistical average of the behavior 
of many, or several, stars, observed by photometry, spec- 
troscopy and asteroseismology. Provided these relations are 
accurate, our predictions are therefore those for a notional 
average field star. Natural scatter means that e.g. not all 
stars hotter or cooler than the Sun will present a smaller 
challenge when it comes to rotation; however, given a suf- 
ficient number of stars, we predict that on average, things 
will be easier. 
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